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Perturbed stellar motions around the rotating black hole in Sgr A* for a generic 

orientation of its spin cixis 

L. loricE 

Empirically determining the averaged variations of the orbital parameters of the stars orbiting 
the Supermassive Black Hole (SBH) hosted by the Galactic center (GC) in Sgr A* is, in principle, 
a valuable tool to test the General Theory of Relativity (GTR), in regimes far stronger than those 
tested so far, and certain key predictions of it like the no-hair theorems. We analytically work 
out the long-term variations of all the six osculating Keplerian orbital elements of a test particle 
orbiting a nonspherical, rotating body with quadrupole moment Q2 and angular momentum S for 
a generic spatial orientation of its spin axis k. This choice is motivated by the fact that, basically, 
we do not know the position in the sky of the spin axis of the SBH in Sgr A* with sufficient 
accuracy. We apply our results to S2, which is the closest star discovered so far having an orbital 
period Ph = 15.98 yr, and to a hypothetical closer star X with Ph = 0.5 yr. Our calculations are 
quite general, not being related to any specific parameterization of k, and can be applied also to 
astrophysical binary systems, stellar planetary systems, and planetary satellite geodesy in which 
different reference frames, generally not aligned with the primary's rotational axis, are routinely 
used. 

PACS numbers: 04.20.-q, 04.70.-s, 04.80.Cc, 91.10.Qm, 91.10.Sp, 97.60. Lf 



I. INTRODUCTION 

There is nowadays wide consensus [l|-y| that the Galac- 
tic Center (GC) hosts a Supermassive Black Hole (SBH) 
d, Q whose position coincides with that of the radio- 
source Sagittarius A* (Sgr A*) [i,0| at d = 8.28 ± 0.44 
kpc from us (8|' for a popular overview of such an object, 
see, e.g., Ref. Q- The Galactic SBH is surrounded by a 
number of recently detected main-sequence stars of spec- 
tral class B [^, [l^l . They have been revealed and tracked 
in the near infrared since 1992 at the 8.2 m Very Large 
Telescope (VLT) on Cerro Paranal, Chile and the 3.58 
m New Technology Telescope (NTT) on La Silla, Chile 
[ll|, and since 1995 at the Keck 10 m telescope on Mauna 
Kea, Hawaii [H. They are dubbed SN, or SO-N in the 
Keck nomenclature, where N is a progressive order num- 
ber. Their relatively fast orbital motions, characterized 
by orbital periods Ph > 16 yr, high eccentricities e > 0.2, 
random orientations i of their orbits in the sky and av- 
erage distances from the SBH r > 2 x lO^r, 
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The direct access to such S/SO stars, and of other 
closer objects which may hopefully be discovered in the 
future, has induced several researchers to investigate var- 
ious predictions that the General Theory of Relativity 
(GTR) directly makes for their orbital motions along 
with other competing effects from standard Newtonian 
gravity which may mask the relativistic ones [14| - |25l |. 



^^g, where rg 



denotes the Schwarzschild radius, allowed to dynamically 
infer a mass of about M « 4 x lO^Mg 0, i, [li for it. 
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Concerning several effects related to propagating elec- 
tromagnetic waves in connection with the stellar orbital 
motions like, e.g., relativistic redshifts, see Ref. |26l - [29| . 
In fact, although the currently known stars, in a strict 
sense, cannot probe the strong field regime of GTR be- 
cause of their relatively large distance from the SBH, on 
the other hand they yield a unique opportunity to put on 
the test GTR in the strongest field regime ever probed so 
far. Indeed, even in the double binary pulsar PSR J0737- 
3039 A/B [30, [31I rg/r is one order of magnitude smaller 
than for S2, which is the closest SBH star discovered so 
far 0. 

In this paper we analytically work out the averaged 
variations of all the six standard osculating Keplerian pa- 
rameters of a test particle caused by the rotation of the 
central object endowed with angular momentum S = Sk 
and quadrupole moment Q2. Note that the stars orbit- 
ing the SBH can safely be considered test particles: their 
masses are about m < 10^^ A/, and relativistic correc- 
tions to their internal structures are assumed to be too 
small to yield noticeable effects on their orbital motions. 
No assumptions about any specific spatial orientation for 
fc are made. Thus, our calculations arc not restricted to 
a particular reference frame, and are valid also for dif- 
ferent scenarios like, e.g., stellar planetary systems and 
planetary satellite geodesy in which natural and/or artifi- 
cial test bodies are employed for testing GTR. Moreover, 
in order to keep our calculations as general as possible, 
we will not adopt any particular representation for k in 
terms of specific angular variables in the sense that we 
will refer it to the global reference frame adopted; for a 
different approach, see Ref. [201 in which k is referred 
to the orbital plane of each star. Concerning the SBH in 
Sgr A* , the orientation of its spin axis is substantially un- 
known, despite the attempts by different groups j32| - |34j 
to constrain it using different parameterizations which 
yielded quite loose bounds. A strategy to partially over- 



come such an obstacle have been recently put forward 
in Ref. [33 ; it is based on the possible observation of 
the apparent shape of the shadow cast by the BH on the 
plane of the sky, and would allow the measurement of S 
and the angle i between k and the line-of-sight. 

The GTR prediction for the standard IPN periastron 
precession, which is analogous to Mercury's well known 
perihelion precession of 42.98 arscec cty~^ and is inde- 
pendent of fc, amounts to about 



w, 



(IPN) 
S2 



45 ± 10 arcsec yr 



(1) 



for S2; the quoted uncertainty comes from the errors in 
the parameters of both the SBH and S2 entering the GTR 
formula: they arc displayed in Table HI The result of Eq. 
([T|). computed in a frame with the SBH at its origin, 
corresponds to a precession ^ as seen from the Earth of 
Cs2 = 27 ± 6 microarcseconds per year (/las yr^^ in the 
following). At present, it is still undetectable from the 
currently available direct astrometric measurements in 
terms of right ascension a and declination S which barely 
cover just one full orbital period of S2. Indeed, according 
to Tabic 1 of Ref. Q , the present-day error in the peri- 
astron is a^ = 0.81 deg = 2916 arcsec over about 16 yr, 
from which an uncertainty of about ct^:, ~ 182 arcsec yr~^ 
in the periastron precession may naively be inferred: it 
corresponds to a limiting accuracy of ct; ~ 110 /xas yr~^ 
in monitoring angular rates as seen from the Earth. As 
wc will see, the sizes of the other precessions of S2 due 
to S and Q2 may be smaller by about 2 and 4—5 or- 
ders of magnitude, respectively for a moderate rotation 
of the SBH. Concerning future perspectives, according 
to Ref. [3g| future astrometric measurements of S2 may 
bring its IPN periastron rate into the measur ability do- 
main; indeed, the periastron advance would indirectly be 
inferred from the corresponding apparent position shift 
in the recorded orbit. Moreover, the ASTromctric and 
phase- Referenced Astronomy (ASTRA) project [33|, to 
be applied to the Keck interferometer, should be able to 
monitor stellar orbits with an accuracy of |38| CAf ~ 30 
/ias as seen from the Earth. The GRAVITY instrument 
[39| , devoted to enhance the capabilities of the VLT inter- 
ferometer (VLTI), aims to reach an accuracy of ctaj ~ 10 
fia.s |39j in measuring astrometric shifts A^ as seen from 
the Earth, which, among other things, would allow ex- 
ploration of the innermost stable circular orbits around 
the SBH 53- 

About testing GTR in the SBH scenario, we make the 
following general considerations. In order to meaning- 
fully compare theoretical predictions for a given effect to 
its empirically determined counterpart, we need to know 
some of the key ambient parameters entering the predic- 
tions independently from the effects themselves we are 
looking for. In the specific case, the mass M, the spin 
S and the quadrupole Q2 of the SBH should be known, 
if possible, independently of the precessions we are going 
to consider. Concerning the SBH mass M, the values 
which we presently have for it can be thought as inferred 
from the third Kepler law used in conjunction with the 



directly measured orbital period Pb, and the semimajor 
axis a empirically determined by modeling the recorded 
stellar orbit in the plane of the sky with an ellipse (see 
Fig. 2 of Ref. Q)- Such a determination of M would be, 
in principle, "imprinted" by GTR itself since it implies a 
correction to the third Kepler law, but it is far too small 
with respect to the present-day accuracy in determining 
Pb- Indeed, it is ap^ ~ 10"^ yr [1, [H, while the IPN 
GTR correction to the Keplerian orbital period is [4l|, \^ 
AP^^™^ ex (37r/c2)VGl/a ~ lO'^ yr for S2. The same 
holds also if M is straightforwardly inferred, in a perhaps 
less transparent manner, as a solvc-for quantity from 
multiparameter global fits of all the stars' data: mod- 
eling or not GTR at IPN level has not yet statistically 
significant infiuence in its estimated values, as shown by 
Tabic 2 of Ref. [g. We stress that, when such an ap- 
proach is followed to test GTR, it is intended that differ- 
ent dynamical models, with and without GTR, arc fitted 
to the same data sets to see if statistically significant 
differences occur in the solve-for estimated parameters. 
The quadrupole moment Q2 of the SBH in Sgr A* may 
be measured, e.g., by means of imaging observations with 
Very Long Baseline Interferometry (VLBI) in the strong 
field regime; see Ref. |43l - |45| for recent reviews and other 
proposals. In regard to the spin S of Sgr A*, one tries 
to gain independent information about S from the inter- 
pretation of some measured Quasi-Periodic Oscillations 
(QPOs) in the X-ray spectrum of the electromagnetic ra- 
diation emitted by the gas orbiting in the accretion disk 
close to its inner edge [4a, |43] . More recent observations 
conducted with the Millimeter Very Long Baseline Inter- 
ferometry (mm- VLBI), probing the immediate vicinity 
of the horizon, have been able to get information on S 
[33, [SJ]- In interpreting such measurements, the validity 
of the Kerr metric [43| as predicted by GTR is assumed, 
thus inferring S from, say, the radius of the inner edge ex- 
tracted from the X-ray diagnostics. It is worth pointing 
out that the mere fact of obtaining a good fit of the Kerr 
metric ^^ to a certain empirically determined quantity 
like, e.g., the X-ray spectrum, getting a reasonable value 
for 5" as a least-square adjustable parameter, cannot be 
considered in itself as a test of the rotation- related predic- 
tions of GTR, also because other competing mechanisms 
to explain, say, the QPOs, whose physics is still rather 
disputed, exist. Independent empirical determinations of 
different effects connected with S are required, and the 
stellar orbital precessions would be just what we need. 
The greater the number of precessions empirically deter- 
mined, the greater the number of GTR tests which can be 
performed. In principle, more than five precessions are 
required since M, S, Q2 and two components of k must 
be determined; see also the discussion in Ref. [20|. Thus, 
the need for more than one star is apparent. Such a num- 
ber of necessary orbital rates may be reduced if some of 
the aforementioned parameters are somehow reliably ob- 
tained from other sources. Of course, also the accuracy 
with which the precessions can be determined plays a 
role, in the sense that the previous reasoning holds in 



the ideal case in which ah the three dynamical effects 
considered are detectable. Basically, it is the same logic 
behind the usual tests in the binary systems hosting at 
least one active radio-pulsar [43 . Indeed, in that case the 
interpretation of just two empirically determined post- 
Keplerian effects in terms of their IPN-GTR predictions 
is not sufficient since it only allows to obtain the masses 
nil and 7712 of the system, which are a priori unknown. 
In the binary pulsar systems the effects which can, actu- 
ally, be inferred from the data are not limited just to the 
post-Keplerian periastron precession. Genuine tests of 
GTR are made when more than two post-Keplerian pa- 
rameters are empirically determined, and the additional 
ones are interpreted with GTR by using in their analyti- 
cal predictions just the previously obtained values for 771 1 
and 7712 figj . 

The plan of the paper is as follows. In Sec. |ll]we review 
basic facts of standard perturbation theory which will be 
applied in Se c. |llI]to Q2 (Sec. UlTl^ and S (Sec. HITbI) . 
In Sec. IIIIDl it is briefly remarked that also gravitational 
waves with ultralow frequency traveling from the outside 
could be constrained by the orbital precessions of the 
stars in Sgr A*. In Sec. IIIICI we compare our results 
to those obtained by Will [20] . Numerical evaluations of 
the effects worked out in Sec. IIIII are presented in Sec. 
IIVI Sec. |V| is devoted to summarizing our findings and 
to the conclusions. 



II. OVERVIEW OF THE METHOD ADOPTED 

Here we deal with a generic perturbing acceleration 
A induced by a given dynamical effect which can be 
considered as small with respect to the main Newtonian 
monopole ^Newton = —GM/r^, where G is the Newto- 
nian constant of gravitation and r is the mutual particle- 
body distance. 

First, A has to be projected onto the radial, trans- 
verse and normal orthogonal unit vectors R, T, N of the 
comoving frame of the test particle orbiting the central 
object. Their components, in Cartesian coordinates of a 
reference frame centered in the primary, are [50| 



cos J? cos u — cos i sin fl sin u 
R = sin n cos u + cos i cos H sin u 
sin i sin u 



— sin u cos f2 — cos i sin i? cos u 
T = \ sin n sin u + cos i cos i? cos u 
sin i cos u 



(2) 



(3) 



the reference {x, y} plane from the reference x direction 
to the intersection of the orbital plane with the reference 
plane {x, y} (the line of the nodes). The angle u ^ uj + f 
is the argument of latitude. In it, uj is the argument 
of pericenter: it is an angle in the orbital plane reck- 
oned from the line of the nodes to the point of closest 
approach, generally known as pericenter. The angle / is 
the true anomaly: lying in the orbital plane, it is counted 
from the pericenter to the instantaneous position of the 
test particle. The angle i is the inclination of the or- 
bital plane to the reference {a;, y} plane. In this specific 
case, we will choose the unit vector p of the line-of-sight, 
pointing from the object to the observer, to be directed 
along the positive z axis, so that the {x, y} plane coin- 
cides with the usual plane of the sky which is tangential 
to the celestial sphere at the position of the BH. With 
such a choice, corresponding to the frame actually used 
in data reduction [SJ, |5l| , i is the inclination of the orbital 
plane to the plane of the sky (i = 90 deg corresponds to 
edge-on orbits, while i = deg implies face-on orbits), 
and n is an angle in it counted from the reference x di- 
rection; it is such a node which is actually determined 
from the observations [^, M, \lM j ^'^'^i i^ general, it is not 
referred to the SBH's equator. Subsequently, the pro- 
jected components of A have to be evaluated onto the 
Keplerian ellipse 



P 



1 + e cos / 



, p^a{l -e^). 



(5) 



where p is the semilatus rectum and a, e are the semima- 
jor axis and the eccentricity, respectively. The Cartesian 
coordinates of the Keplerian motion in space are [50| 



' (cos f] cos u — cos i sin J? sin u) . 
■ (sin J? cos u + cos i cos i? sin u) , 



(6) 



rsim smu. 



Then, Aji^At^An are to be plugged into the right- 
hand-sides of the Gauss equations for the variations of 
the osculating Keplerian orbital elements [42, ISJ] . 

Their right-hand-sides, computed for the perturbing 
accelerations of the dynamical effect considered, have to 
be inserted into the analytic expression of the time vari- 
ation d W /dt of the osculating Keplerian orbital element 
^ of interest. Then, it must be averaged over one orbital 
revolution by means of [52| 



dt 



(1 - e^)3/2 
n(l + ecos/) 



r#, 



(7) 



TV: 



sin i sin J? 

— sin i cos J? 

cosi 



(4) 



In Eq. (|3])-Eq. ^, the angles fi,u,i are as follows. The 
angle /? is the longitude of the ascending node: it lies in 



where n = ^/GM/a? is the Keplerian mean motion re- 
lated to the orbital period hy n — 27r/Fb, to obtain 
{dW/dt). As a general remark, we point out that it 
would be incorrect to make inferences about the aver- 
aged orbital effects (d^/dt) from a simple inspection of 
the analytic form of the components Aj^^ At^ An of a 



given perturbing acceleration A, apart from simple triv- 
ial cases. The actual calculation must ultimately be per- 
formed in full as previously outlined. Indeed, it may well 
happen that nonzero components of A yield vanishing 
averaged variations (d'P/dt) for some Keplerian orbital 
elements (?". A trivial case occurs, of course, when one or 
more components of A are identically zero. Conversely, it 
would be incorrect to argue that certain components of A 
should necessarily vanish only because the averaged vari- 
ations of the orbital elements involving them have been 
found to be zero. Moreover, simple back-of-the-envelope 
numerical evaluations of the size of the averaged vari- 
ation {d'F/dt) of a given Keplerian orbital element 'P 
which are based on the order of magnitude of A may be 
misleading as well. Indeed, it may happen that the final 
result (dW/dt) of the complete calculation retains a mul- 
tiplicative factor e", a — ±1, ±2, ±3, . . . which can cause 
a notable quantitative difference with respect to what 
naively guessed, especially for low-eccentricity systems. 



III. CALCULATION OF THE LONG-TERM 
ORBITAL EFFECTS 



J53l 1541 ] , an electrically neutral BH is completely charac- 
terized by its mass M and angular momentum S only. 
As a consequence, all the multipolc moments of its ex- 
ternal spacetime are functions of M and S [5^, [Sg] . In 
particular, the quadrupole moment of the BH is 



^2 = 






(9) 



The spatial orientation of the BH's spin axis can be con- 
sidered as unknown. Thus, looking for a more direct 
connection with actually measurable quantities, we will 
retain a generic orientation for k in the ongoing calcula- 
tion, i.e., we will not align it to any of axes of the refer- 
ence frame used. After having computed the R — T — N 
components of Eq. ^ by means of Eq. ©-Eq. (g]) as 



R 

(Q2 

T 

(Q2 

N 



A(Q2) . TV, 



(10) 



A. The long-term precessions caused by the 

quadrupole mass moment of the central body for an 

arbitrary orientation of its spin axis 

The acceleration experienced by a test particle orbit- 
ing a nonsjjherical central mass rotating about a generic 
direction k is 



^(Q2) ^ _2^ 



302 

2r4 



1-5 f -fc 



2{f-k\k 



(8) 

where Q2 is the quadrupole moment of the body, with 
[Q2] = L^T^^. A dimensionless quadrupole parameter 
J2 can be introduced by posing Q2 -> GA/T^.^ J2, where 
Ti-e is the equatorial radius of the rotating body. Ac- 
cording to the "no-hair" or uniqueness theorems of GTR 

I 



to be evaluated onto the unperturbed Keplerian ellipse, 
it is possible to obtain 



(11) 



for the semimajor axis and the eccentricity, as in the 
standard calculations [52| in which k is usually aligned 
with the z axis. 

Instead, the inclination i undergoes a long-term varia- 
tion given by 



3g2 



2VGM^{l-e'^y 



1 [n,i\k] , 



(12) 



with 



3 I J7, i; fc j = \kx cos H + ky sin J? 



k^, cos i + sin i ( fc.^ sin il —- ky cos il 



(13) 



r 



If kx ^ ky = 0, as in the usual calculation [52| . i stays Concerning the node Q, its long-term variation is 

constant. 

'dn\ SQ2 



dt / AVGMa^il~e^y 



(/2,i;fc), (14) 



with 



D ( n,i;k] = 2kz cos 2i esc i ( kx sin fi — ky cos J? ) + cos i 



(kx su 



f^x + ^y 2k~ 



[ky k^j 



cos 2/2 — 2kxky sin 2i7 



(15) 



Notice that kx = ky = in Eq. (fTS]) yields the standard secular precession [521 with 

D (i) = — 2cosz. 



(16) 



The long-term change of the argument of pcricentcr uj with 
is a httle more cumbersome. It is 



duj 



3Q2 



leVGA/a^ (1 - e^f 



( i7,i; fc 



(17) 



{n,i;k) 



= 8 - llkl - llkl - 2kl + (fc2 - kl ) cos 2/2- 



2kz (cot i — 5 cos 3i esc i) I ky cos il ^ k^ sin J? 



5 cos 2i 



2kxky sin2j?- 



2fc2 - fc2 - A:2 -h f A:^ - A:^ j cos 21? -f- 2k^ky sin 21? 



(18) 



In the case fcr 



Eq. HH]) reduces to 



by 



o{i) =:2(3 + 5cos2i) =4(4-5sin^ 



(19) 



which yields the standard expression for the secular pre- 
cession of the pericenter [52| . 

The longitude of the pericentre ru = oj + H, which is a 
"dogleg" angle, experiences a long-term variation given 

I 



/ dzu 



3Q2 



leVCMa^ {1 - e^y 



^ (n,i;k 



(20) 



with 



^(n,i;k] = 8-nkl-nkl~2k'^+ (kl + kl -2kfj {4 cosi- 5 cos 2i)- 

- 4:(kl-kl\ (3-^5cosi)sin2(i)cos2/2-2fcyfc^sec(f) [sin(f) -|- 5sin (f )] cos /?+ (21) 

-I- 2k^kzsec{^) [sin(^) -h5sin(f )] sin ~ Sk^ky sin^ (|) (3 + 5 cos i) sin 2/2. 

I 



For kx = ky = Eq. ([2T|) reduces to is 

QJ(i) = 2[3- (4cosi-5cos2i)] = 4(4-5sin^ i-2cosi), 

(22) 
which yields the usual expression for the secular rate of 

TO dn. 

Finally, the long-term change of the mean anomaly 
Ai = n{t — ip), where ip is time of passage at pericenter, with 

I 



/dM\ 
\dt / 



3Q2 



WjGMa'^ (1 - e^y 



m[n,i;k], (23) 



m 



(f2,i;k\ = -8 + 9fc2 + 9fc2 + Qk^ + 3 (kl + kl - 2kf\ cos 2i + 6 (k^ - fc^^j sin^ i cos 2/2- 



12 



( 



kz sin 2i ( ky cos fi — k^ sin /? J -I- fca;A;T, sin i sin 2 J? 



) 



(24) 



Also in this case, for k^ = ky = the standard secular of Eq. (fTT|) -Ea. p4)) is not limited just to the BH arena, 

being then generally valid also for astrophysical binary 
systems, stellar planetary systems, and planetary satel- 
lite geodesy. In particular, they could be useful when 



precession [52| is recovered since Eq. ([24|) reduces to 
m {i) = -2 (1 -f 3 cos 2i) = -4 (2 - 3 sin^ i) . (25) 
Incidentally, we remark that the field of applicability 



satellite-based tests of GTR are performed or designed 
(See See. V). 



B. The Lense-Thirring long-term precessions for a 

generic orientation of the spin axis of the central 

body 

According to GTR, the gravitomagnetic Lense- 
Thirring acceleration felt by a test particle moving with 
velocity v around a rotating body with angular momen- 
tum S = Sk at large distance from it is 



(LT) 



(- 

V c 



xB„ 



(26) 



In Eq. p6|) the gravitomagnetic field Bg, far from the 
central object where the Kerr metric [4^ reduces to the 
Lense-Thirring one, is 



B„ 



GS 



f,j.O 



fe — 3 I fe • r I r 



(27) 



Concerning S*, the existence of the horizon in the Kerr 
metric |48| implies a maximum value for the angular mo- 
mentum of a spinning BH [57|, |58[, so that S = XgS'max, 
with 



•^max 



APG 



(28) 



If Xg > Ij the Kerr metric [4g would have a naked sin- 
gularity without a horizon. Thus, closed timelike curves 



could be considered, implying a causality violation J59| . 
Although not yet proven, the cosmic censorship con- 
jecture [60| asserts that naked singularities cannot be 
formed via the gravitational collapse of a body. If the 
limit of Eq. ((28| is actually reached or not by astrophys- 
ical BHs depends on their accretion history |6l| . In the 
case of Sgr A*, it may be Xg ~ 0.44-0.52 [46l.l47| or even 
less [3j, [SJ]. Contrary to BHs, no theoretical constraints 
on the value of Xg exist for stars. For main-sequence 
stars, Xg depends sensitively on the stellar mass, and can 
be much larger than unity [63 - l65l | . The case of compact 
stars was recently treated in Ref. [6g|, showing that for 
neutron stars with M > IMq it should be Xg ^ 0.7, in- 
dependently of the Equation Of State (EOS) governing 
the stellar matter. Hypothetical quark stars may have 
Xg > 1, strongly depending on the EOS and the stellar 
mass l66i. 



In the standard derivations of the Lense-Thirring effect 
pl\ existing in literature the reference {x, y} plane was 
usually chosen coincident with the equatorial plane of the 
rotating mass. In principle, the Lense-Thirring orbital 
precessions for a generic orientation of S could be worked 
out with the Gauss equations in the same way as done 
for Q2. Anyway, they were recently worked out [68[, in a 
different framework, with the less cumbersome Lagrange 
planetary equations [52] ■ For the reader's convenience, 
we display here the final result 



/ da\ 
\dt / 

(f) 

(i) 

/ dn\ 

\ dt / 

/ dui\ 
\ dt / 



/rfro\ 
\ dt / 

/ dM\ 
\ dt / 



0, 
0, 

2GS(fcrr COS n + ky sin fi) 
c2a3(l-e2)3/2 ' 

2GS\kz +cot i{ky cos Q — kj^ sin 12 j 

c2a3(l-e2)3/2 ' 

GS\Qkz cos z+(3 cos 2i—l) esc i{ky cos H — k^ sin H) 

c2a3(l-e2)3/2 ' 

GS j 4 fc- cos z-t-sin iik^ sin O — ky cos f2) — 2 /c- sin i-t-cos iiky cos O — kj. sin Oj tan(2/2) j- 



(29) 



c2Q3(l_e2)3/2 



0. 



Notice that Eq. (P^ yields just the usual Lense-Thirring 
secular rates [42,|63] for kx = ky = 0. Contrary to such a 
scenario, the inclination i experiences a long-term grav- 
itomagnetic change for an arbitrary orientation of S: it 
is independent of the inclination i itself. 



C. A comparison with a different approach 



Will (20[ refers k to the orbital plane of a generic star 
by choosing ep,eq^h as orthonormal vectors: Sp is di- 
rected along the line of the nodes, e^ lies in the orbital 



plane perpendicularly to Ep, and h is directed along the 
orbital angular momentum. The unit vectors e^ and Cg 
can be obtained from Eq. ([2]) and Eq. (j3|), respectively, 
by posing u -^ 0, while h coincides with Eq. (U). Thus, 
one has 



(fch 



kx = kp cos n + ikh sin i — kq cos i J sin i7. 



fcj, = fcp sin i? — I fc/i sini — fcg cosi 1 cos i?, ion\ 

kz ~ kh cos i + kq sin i. 



D. Stellar orbital perturbations caused by ultralow 
frequency gravitational waves 

The stars orbiting the SBH in Sgr A* could also 
be used, in principle, as probes for detecting or con- 
straining plane gravitational waves of ultralow frequency 
(j/ PS 10~® — 10~^° Hz or less) impinging on the sys- 
tem from the outside. Indeed, the passage of such waves 
through the orbits of the closest stars would cause long- 
term variations of all their Keplerian orbital elements, 
apart from the semima .ior axis a. They have recently 
been worked out in Ref. [6^ for general orbital configura- 
tions, i.e., without making a-priori assumptions on their 
inclinations and eccentricities of the perturbed test par- 
ticle, and arbitrary directions of incidence for the wave. 
Conversely, gravitational waves can be generated within 
the stellar system of Sgr A* , as discussed in Ref. [70] ■ 



IV. NUMERICAL EVALUATIONS 



where kp,kq,kfi can straightforwardly be expressed in 
terms of the polar angles a and (3 used by Will [201 ^^ 
the orbital frame. 



Inserting Eq. (pO| into the equations of Sec. IIII Al and 
Sec. IIII B I allows us to obtain Eq. (2a), Eq. (2b), and 
Eq. (2c) of Ref. [2^ after some algebra. 



In Tabic |T] we quote the relevant physical and orbital 
parameters for the SBII-S2 system. The orbital period 
of S2 is Pb = 15.98 yr, so that the astromctric measure- 
ments currently available cover a full revolution of it. 

The quadrupole-induced precessions of Eq. (fT2|) -Eq. 
p4| are all linear combinations of the products of the 
components of k plus, sometimes, a term independent of 
fc: they can be cast into the form 



~dt 



Dq (Q2, a, e, J, ri) + - ^ Dsi (Q2, a, e, i, il) kgh, s,l = x,y,z, ^ ^ i, i},u>,M. 



(31) 



The numerical values of the coefBcients Dq and Dsi = 
Dig for S2, in /xas yr~^, are quoted in Table |lTl The 
largest effects occur for oj and M. because of Dq, which 
is of the order of « 1 milliarcsec yr~^ (mas yr~^). The 
other terms are damped by the square of the components 
of fc. Moreover, partial mutual cancellation may occur 
depending on the orientation of the SBH spin axis. 

The Lense-Thirring precessions of Eq. pQ]) are all lin- 
ear combinations of the components of fe: they can be 
cast into the form 



'i)-T.c.is. 



a,e,i,fl)kj, j = x,y,z, 'l' = i,Q,,uj. 

(32) 
The numerical values of the coefficients Cj for S2, in arc- 
^ , are listed in Table IIIII They are of the order of 
i.e., orders of magnitude larger 
than the quadrupole precessions of Table |llj also in this 
case, partial mutual cancellations may occur depending 
on fc, thus impacting the detectability of the gravitomag- 
netic rates. 



sec yr 

about 10^^ arcsec yr~^ 



I 

The figures of Table [H] and Table IIIII can be compared 
with the present-day accuracies in empirically determin- 
ing the orbital precessions of S2 listed in Table IIVI They 
are of the order of 10^ — 10^ arcsec yr~^. The Lense- 
Thirring precessions of S2 (Table IIIip are about three 
orders of magnitude smaller than the current accuracy, 
while the quadrupole effects of Table |ll] are negligibly 
small. 

By considering a fictitious star X with, say, the same 
orbital parameters of S2 apart from the semimajor axis 
a, assumed to be one order of magnitude smaller so that 
its orbital period would just be Pb = 0.5 yr, it turns out 
that its IPN GTR periastron precession would be as large 
as 4 deg yr~^, while its Lense-Thirring and quadrupole 
precessions would be of the order of about « 10^ arcsec 
yr~^ and w 1 arcsec yr"^, respectively. 

If, as expected, angular shifts of A^ « 10 ^as, as 
seen from the Earth, will really become measurable in 
future thanks to GRAVITY and ASTRA, this would im- 
ply an accuracy of the order of A <?" w (d/a) A^ — 16 
arcsec for S2, and 160 arsec for a star one order of mag- 



TABLE I. Relevant physical and orbital parameters of the SBH-S2 system in Sgr A* (first row), and their uncertainties (second 
row) . The Keplerian orbital elements of S2 were retrieved from Table 1 of [g| . The figure for Xg comes from Ref. [43| , while the 
one for the gravitational parameter /i = GM is from a multi-star fit yielding /i = (4.30 ±0.50) x 10^/^0 [g]. The quoted value in 
m for the semimajor axis of S2 was obtained by multiplying its angular value a — 0.1246 ± 0.0019 arcsec Q] by the distance of 
the SBH d — 8.28 ±0.44 kpc [g|. For the angular momentum and the quadrupole moment of the SBH we used S — XgiM^G^jc 
and Qi = -{S'^G)l(i^M) = -xl{G^ M^) / c^ . The orbital period of S2 is Pb = 15.98 yr = 5.04 x 10* s. The figures for S 
and Q2, obtained in the hypothesis that GTR is correct, strongly depend on Xs, which is, at present, highly uncertain. For 



example, the authors of Ref. 



yield Xg = 0.44 ±0.08, while for the authors of Ref . [331. |34| it could be even smaller. We will 



use them to indicatively give order-of-magnitude evaluations of the additional orbital precessions which would occur because 
of 5* and Q2 according to GTR. 



^ (m^ s-^) 



S (kg m^ s-^) 



(m" 



Xg 



a (m) 



(dcg) 



n (deg) 



5.70 X 10^ 
6.6 X 10^^ 



8.46 X 10" 
4.66 X 10^ 



-6.22 X lO'* 
6.58 X 10*^ 



0.52 
0.26 



1.54 X 10" 
8 X 10^^ 



0.8831 
0.0034 



134.87 
0.78 



226.53 
0.72 



TABLE IL Coefficients of the quadrupole precessions of S2, 
in /^as yr~^, according to Table|I] GTR was assumed for Q2, 
with Xg ~ 0.52. 





Do 


D^2 


D,. 


D,2 


D^y 


D^, 


Dy. 


i 





406 


-406 





43 


558 


588 


n 





427 


384 


-810 


-809 


-5 


5 


UJ 


-1149 


1568 


1584 


294 


293 


-1254 


1189 


M 


-539 


595 


616 


406 


405 


-587 


556 



TABLE III. Coefficients of the Lense-Thirring precessions of 
S2, in arcsec yr~^, according to Tabled In particular, Xg ~ 
0.52 was used for the spin of the SBH. 



Ox 



Oil 



a- 



UJ 



-0.14 
-0.15 
0.11 



-0.15 

0.14 

-0.10 





0.21 

0.45 



nitude closer to the SBH. If such targets will be dis- 
covered, their Lense-Thirring shifts should become de- 
tectable after some years, while the (52— induced per- 
turbations would still remain hard to measure, even for 
ewO.9. 



TABLE IV. Naive evaluations of the uncertainties in the secu- 
lar variations of the S2 osculating Keplerian orbital elements, 
in arcsec yr^^, obtained by dividing the errors in the ele- 
ments from Table 1 of Ref. [g| by a time interval AT ~ Pb- 
Concerning the mean anomaly, its uncertainty was evaluated 
from that of the time of periastron passage tp, released in 
Ref. [g , according to the expression for the mean anomaly at 
the epoch of periastron passage Mo = —ntp; also the errors 
coming from a and jj, through n were taken into account. 



176 



163 



182 
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SUMMARY AND CONCLUSIONS 



We analytically worked out the long-term, i.e., aver- 
aged over one full revolution, variations of all the six oscu- 
lating Keplerian orbital elements of a test particle orbit- 
ing a nonspherical, spinning body endowed with angular 
momentum S and quadrupole moment Q2 for a generic 
spatial orientation of its spin axis k. We did not restrict 
ourselves to any specific orbital configuration of the par- 
ticle. Here we applied our results to the stars orbiting 
the SBH in Sgr A*: those identified so far are moving 
along highly elliptical trajectories with periods Pb > 16 
yr. The current level of accuracy in empirically deter- 
mining the precessions of the angular orbital elements of 
S2, having Pb = 16 yr, can be evaluated to be of the or- 



der of « 10^ — 10"^ arcsec 



yr 



The predicted IPN GTR 



periastron precession of S2, which is independent of the 
orientation of the spin axis of the SBH, is 40 ± 10 arcsec 
The predicted GTR spin and quadrupole-induced 



yr" 



^ arcsec yr ^ 



precessions of S2 are of the order of ~ 10 
and sa lO'^ — lO'^/ias yr~^, respectively: they depend on fc, 
and partial cancellations among their components may 
occur, thus reducing their magnitude. Concerning hy- 
pothetical stars with orbital periods of less than 1 yr, 
not yet discovered, the IPN GTR periastron precessions 
would be as large as some deg yr~^, while the 5* and Q2 
effects would be of the order of w 10^ arcsec yr~^ and 
« 1 arcsec yr~^, respectively. Planned improvements of 
the infrared telescopes used so far aim to reach an accu- 
racy level of « 10 /ias at best in measuring angular shifts 
as seen from the Earth corresponding to stellar orbital 
shifts of about 1.6 x 10^ — 10^ arcsec for S2 and stars 
closer than it by one order of magnitude, respectively. 
Finally, we stress that our calculations are not restricted 
to any specific coordinate system. Thus, in the form in 
which we obtained them, they can fruitfully be used also 
in other scenarios like astrophysical binaries, stellar plan- 
etary systems and planetary satellite geodesy in which 
tests of post-Newtonian gravity may be involved. Exam- 
ples of that can be found, e.g., in Ref. [7l[, pertaining 
the Lense-Thirring tests with the LAGEOS satellites in 
the terrestrial gravitational field, and in Ref. [T^l dealing 



with Mercury and the gravitomagnetic field of the Sun. 
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